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1.0  INTRODUCTION 


The  purpose  of  this  report  is  to  document  a  radar  sea  clutter  model  for  inclusion  as  a  Navy  standard 
model  in  the  Oceanographic  and  Atmospheric  Master  Library  (OAML). 

The  generalized  form  of  the  radar  equation  can  be  written  as 


Pr 


p,gWl} 

(4  nyR* 


oeF 4 


(W) 


(1) 


where  Pr  is  power  received,  Pt  is  transmitted  power,  G  is  antenna  gain,  X.  is  radar  wavelength,  Ls  is  system 
losses.  R  is  range,  Oc  is  average  clutter  cross  section,  and  F  is  the  pattern  propagation  factor.  If  the 
antenna  pattern  effects  are  separated  from  F,  Eq.  1  becomes 
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where /(a)  is  the  antenna  pattern  factor  for  angle  a,  and  F  is  now  the  propagauon  factor.  Sea  clutter  cross 
section  can  be  factored  further  by 


Oc  =  (P  Ac 


(m2) 


(3) 


where  o°  is  normalized  radar  reflectivity  of  the  sea  surface  (m2/m2),  and  Ac  is  the  radar's  clutter 
resolution  cell 
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cr 

2 


9 

2  ln(2) 


R  sec(v) 


(m2) 


(4) 


In  Eq.  4,  c  is  the  speed  of  light,  r  is  the  radar's  compressed  pulse  length  (or  pulse  length  if  pulse 
compression  is  not  used),  6  is  the  horizontal  beamwidth  (one-way,  3  dB),  and  is  grazing  angle  at  the 
surface.  For  angles  of  interest  here,  V  ^  10  degrees  and  sec  (tP)  *■  1.  In  the  type  of  semiempirical  sea 
clutter  model  being  proposed  here,  the  propagation  factor  is  considered  to  be  included  in  the  models  for 
oc  ;  that  is,  cPacF*  can  be  thought  of  as  effective  cross  section,  °c.  Thus,  the  radar  equation  is  now 


PiG1X1Lsf(a) 

(4  ji)3R*  °c 


A 


(W) 


(5) 


and  all  of  the  environmental  and  propagation  considerations  are  contained  in  the  model  of  o°  .  An  existing 
normalized  sea  clutter  model'  has  been  shown  capable  of  approximating  observed  clutter  data  if  grazing 
angles  are  calculated  based  upon  typical  surface-layer  refractive  conditions  over  the  ocean.2  The  model 
proposed  here  combines  this  adaptation  of  the  (/’model  with  the  radar  resolution  cell,  Ac  ,  to  provide  a 
model  of  oc  in  Eq.  5  versus  range  under  varying  propagation  conditions  for  typical  radar  frequencies. 

In  section  2.0,  the  Inputs,  outputs,  and  limits  of  the  sea  clutter  model  are  defined.  In  section  3  0,  the 
model  is  described  in  detail.  In  section  4.0,  test  data  are  shown  to  verify  the  proper  operation  of  the 
computer  program  listed  in  the  appendix. 
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2.0  INPUTS,  OUTPUTS,  AND  LIMITS 


Environmental  measurement.*:  and  radar  system  parameters  are  required  to  calculate  the  average  sea 
clutter  radar  cross  section  as  the  output.  The  environmental  input  parameters  are  summarized  in  Table  1. 


Table  1.  Environmental  input  parameters. 


Parameter 

Limits 

Precision 

Evaporation  duct  height 

0  to  40  m 

0.1  m 

Modified  refracuvity  profile 

0  to  10,000  m 

1  m 

0  to  2000  M 

1  M 

Wind  speed 

0  to  25  m/s 

1  m/s 

Wind  direction  relative  to 

0  to  180  deg 

5  deg 

antenna  boresight 

2.1  ENVIRONMENTAL  PARAMETERS 

2.1.1  Evaporation  Duct  Height 

The  evaporation  duct  height  can  be  calculated  from  measurements  of  wind  speed,  air  temperature, 
relative  humidity,  and  sea  surface  temperature,3  or  derived  from  a  statistical  data  base.4  The  evaporation 
duct  height  is  used  to  generate  a  modified  refractivity  profile  for  the  atmospheric  surface  layer.  The 
M-profile  is  subsequently  used  in  determining  grazing  angles  at  the  sea  surface  for  frequencies  of  2  GHz 
and  greater.  The  evaporation  duct  height  is  limited  to  0  to  40  meters. 

2.1.2  Modified  Refractivity  Profile 

The  modified  refractivity  profile  is  represented  by  couplets  of  height  in  meters  and  refractivity  in 
M-units  from  the  surface  to  some  height  greater  than  the  radar  altitude.  This  profile,  which  can  be  derived 
from  radiosonde  or  refractometer  measurements  or  a  statistical  data  base,4  is  used  in  the  model  to 
calculate  an  effective  earth  radius  for  the  determination  of  grazing  angles  at  the  sea  surface  for  frequencies 
less  than  2  GHz. 

2.1.3  Wind  Speed  and  Direction 

Wind  speed  and  direction  can  be  measured  by  any  typical  cup  or  aerovane  anemometer.  The 
measurement  should  be  made  in  a  location  that  is  representative  of  the  wind  flow  over  the  ocean  at  the 
standard  anemometer  height  (19.5  meters).  The  valid  range  of  the  windspeed  is  0  to  25  m/s.  Wind 
direction  for  the  sea  clutter  model  is  defined  relative  to  the  direction  in  which  the  radar  antenna  is 
pointing;  the  measured  wind  direction  must  be  converted  to  this  definition  in  the  range  of  0  to  180 
degrees.  Wind  speed  is  used  to  calculate  a  wind  speed  factor  in  the  sea  clutter  model.  Under  the 
assumption  of  a  fully  risen  sea  and  steady-state  conditions,  wind  speed  and  wave  height  are  highly 
correlated,  and  wind  speed  is  used  to  calculate  an  average  wave  height  for  the  interference  factor.  Under 
transient  conditions,  wave  height  and  wind  speed  may  not  be  well  correlated,  and  the  model  could  be 
modified  to  accept  both  wind  speed  and  average  wave  height  as  inputs. 
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Wind  direction  relative  to  the  radar  antenna  boresight  is  used  to  calculate  an  upwind/downwind  factor 
that  varies  from  a  maximum  looking  into  the  wind  to  a  minimum  looking  with  the  wind.  An  intermediate 
clutter  level  calculated  in  the  crosswind  direction  (90  degrees)  yields  a  level  that  is  *■  ±  5  dB  of  the  clutter 
level  in  the  upwind/downwind  directions. 

2.2  RADAR  PARAMETERS 

The  required  radar  system  parameters  are  summarized  in  Table  2.  Antenna  height,  polarisation,  and 
radar  frequency  are  used  in  the  calculations  of  normalized  radar  cross  section,  <J°.  The  horizontal 
beamwidth  and  compressed  pulse  length  are  used  in  the  calculation  of  the  radar  clutter  resolution  cell.  If 
the  radar  receiver  does  not  use  pulse  compression,  the  pulse  length  is  used. 

Table  2.  Radar  system  parameters. 


Parameter 

Limits 

Precision 

Antenna  height 

1  to  10000  m 

1  m 

Radar  frequency 

100  to  20000  MHz 

1  MHz 

Antenna  polarization 

H,  V,  C 

- 

Horizontal  beamwidth 

>0  to  45  deg 

0.1  deg 

Compressed  pulse  length 

>0  to  1000  *is 

0.1  jis 

Maximum  range 

500  km 

1  km 

*  H  =  horizontal,  V  =  v..  jcal,  C  =  circular. 


2.3  OUTPUT 

The  output  of  the  sea  clutter  model  is  clutter  level  in  decibels  relative  to  a  1  m2  target  versus  range  in 
kilometers. 
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3.0  MODEL 


The  sea  clutter  model  is  composed  of  two  main  functions.  First,  grazing  angle  at  the  sea  surface,  .  is 
determined  from  inputs  of  evaporation  duct  height,  6  ,  effective  earth  radius,  ac  ,  radar  antenna  height, 
Ht ,  and  radar  frequency,  /.  Next,  normalized  radar  cross  section,  a°,  is  calculated  from  inputs  of  grazing 
angle,  wind  speed,  Vw ,  wind  direction  relative  to  the  antenna.  0 ,  radar  frequency,  and  antenna 
polarization.  The  normalized  radar  cross  section  can  then  be  used  directly  to  calculate  average  radar  cross 
section,  ac ,  versus  range  out  to  the  maximum  range  specified. 


3.1  GRAZING  ANGLE  AT  SEA  SURFACE 


At  radar  frequencies  less  than  2000  MHz,  the  effective  earth  radius  is  assumed  to  be  the  dominant 
factor  in  determining  grazing  angle  versus  range.  The  minimum  grazing  angle  allowed  is  1.745 
X  10‘3  radians  (0.1  degree);  the  maximum  grazing  angle  is  1.745  x  10*'  radians  (10  degrees).  The 
range,  r,  in  kilometers  at  which  a  specific  grazing  angle  occurs  is  determined  from  the  quadratic  formula 


r 


-t>  +  Jb2-4ac 
la 


(km) 


(6) 


with 


a  =  l 

b  =  la,\f> 
c  =  -2a,//,/ 1000 


where  a,  is  in  kilometers,  tp  is  in  radians,  and  //,  is  in  meters.  The  effective  earth  radius  is  determined 
by  using  subroutines  MFROF,  INSRT,  PUSH,  and  GETK  from  the  standard  propagation  model.6  The 
elevation  angle,  or,  at  the  antenna  that  corresponds  to  a  given  range  is 


H,  _  J_ 
lOOOr  la. 


(rad) 


(7) 


Grazing  angles  at  ranges  intermediate  to  those  corresponding  to  the  maximum  and  minimum  grazing 
angles  are  determined  by 


JL L.  +  JL 

1000r  2a, 


(rad) 


(&) 


where  range  increments  are  1/100  of  the  maximum  range. 

At  radar  frequencies  of  2000  MHz  and  greater,  the  evaporation  duct  can  have  a  significant  effect  on 
the  variation  of  grazing  angle  versus  range.  Instead  of  an  effective  earth  radius  being  used  in  the  ray*optic$ 
calculations,  a  refractive  profile  for  the  evaporation  duct  must  be  generated.  The  M-v alue  at  the  sea 
surface,  M0,  is  arbitrarily  set  equal  to  340,  and  the  profile  is  calculated  from 


Af(z)  =  M0  +  0.125?  -  0.125dln 


(9) 
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where  In  is  the  natural  logarithm  and 
2  =  «*.*  =  -2,-  1,0,  1.2. 3.  4 

for  heights,  r,  above  the  surface.  The  aerodynamic  roughness  of  the  sea  surface.  20.  is  assumed  constant  at 
0.00015  meter.  The  minimum  value  in  the  evaporation  duct  profile  occurs  at  the  evaporation  duct  height 
(where  2-6)  and  must  also  be  included  in  the  profile.  (Note  that  Eq.  9  can  be  used  to  generate  a  profile 
with  subre tractive  layers  if  a  negative  <5  is  entered.  The  physical  meaning  of  negative  6  is  the  height  at 
which  dM/dz  reaches  a  value  of  +  0.250M/m.  However,  the  evaporation  duct  model3  has  been  arbitrarily 
limited  to  positive  duct  heights.)  To  complete  the  profile  and  as  a  convenience  for  ray-optics  calculations, 
the  height  of  the  radar  is  included  in  the  height  array  and  the  Af-v alue  at  the  height  of  the  radar  is 
included  in  the  M-array.  This  Af-value  is  determined  by  linear  interpolation  between  adjacent  levels  or 
extrapolation  at  0.118  M/m  above  the  last  level  in  the  profile. 

The  minimum  elevation  angle  allowed  is  -1.745  x  10" 1  radians  (-10  degrees).  The  maximum 
elevation  angle  is  calculated  from 

=  -[2  x  10-‘(MA(-A/mln)]’/i  -  1  x  10-‘  (rad)  (10) 

where  M^t  is  the  value  of  M  at  the  radar  antenna  height  and  Mmin  is  the  minimum  value  of  Afm|n  in  the 

profile  below  the  height  of  the  radar.  If  the  radar  height  is  less  than  the  evaporation  duct  height,  then 

a  m«  =  -C2  x  10-‘(M*,  -  Af,)]1'2  -  1  x  lO'6  (rad)  (11) 

where  Mt  is  the  value  of  M  at  the  evaporation  duct  height.  The  grazing  angle  corresponding  to  each 
elevation  angle  is 

y>  =  [pt1  -  2  x  10 -\Mh,  -  A/0)],/2  (rad)  (12) 

The  elevation  angle  is  gradually  decreased  between  the  maximum  and  minimum  angles  according  to 

A  a  =  a.miD  (rad)  (13) 

d 

where  d  has  a  value  dependent  upon  the  last  value  of  \p  : 

y  <  0.0175  d  -  200 

0.0175  S  ^<0.08  d=  20 

d  -  4 

V»  S  0.08 

The  range  corresponding  to  each  elevation  angle  (and  grazing  angle)  is  determined  by  tracing  a  ray  from 
the  radar  antenna  height  through  the  successive  layers  of  the  Af-profile  to  the  surface.  For  downgoing  rays 
(a<0).  the  angle,  or’,  at  the  base  of  the  layer  is 

a  '  =  -  [a2  +  2  x  10-6dAta/i(/tM  -  fj,)]1/2  (rad)  (14) 
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and  the  range,  r',  traveled  is 


r' 


a'  -  a 

r  +  1  x  lO~}dMdh 


(rad) 


05) 


where  dMdh  is  the  refractive  gradient  (M/m)  of  the  layer.  If  the  radar  antenna  height  is  less  than  the 
evaporation  duct  height,  the  rays  are  initially  upgoing  and  are  traced  to  the  level  of  the  evaporation  duct 
height.  For  upgoing  rays,  the  elevation  angle  at  the  top  of  the  layer  is 


a'  ~  [a1 +  2  x  lO-6dMdh(hltl-h,)]in 


(rad) 


(16) 


and  r'  is  calculated  as  above.  Once  the  range  traveled  by  the  ray  to  the  height  of  the  evaporation  duct  is 
known,  it  is  doubled  and  ray  tracing  continues  with  downgoing  rays  at  the  radar  antenna  height  (Fig.  1). 


Figure  1.  Ray  trace  for  antenna  (a)  above 
the  evaporation  duct  height  and  (b)  below 
the  evaporation  duct  height. 

These  calculations  of  the  grazing  angle  do  not  take  into  account  the  effect  that  a  surface-based  duct 
from  elevated  refractive  layers  could  have  on  the  clutter  return  for  a  radar  in  the  duct. 

3.2  RADAR  CLUTTER  CROSS  SECTION 

Wind  speed  and  direction  and  average  wave  height  are  used  to  describe  the  condition  of  the  sea 
surface.  Radar  frequency  and  polarization  are  system  parameters  that,  along  with  the  grazing  angle, 
describe  the  effects  of  the  interaction  between  an  electromagnetic  wave  and  the  sea  surface.  Radar 
wavelength,  A,  is 
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x  122. 

/ 


(m) 


(17) 


where  /  is  frequency  in  megahertz.  Average  wavelength  of  the  ocean  windwaves,  /ta»,  is 


where  V*  is  wind  speed  in  m/s.  Equation  18  assumes  a  fully  risen  sea  and  steady-state  conditions.  Under 
transient  conditions,  Vw  and  ha,  are  decoupled;  the  model  can  be  used  with  these  two  variables  input 
separately.  Three  factors  need  to  be  calculated.  The  interference  factor  is 


A/  = 


1  +oJ 


(19) 


where 


(14.4  X  +  S  Styhai 

- I - 

The  upwind/downwind  factor  is 

Au  =  exp|O.2cos(0)(l-2.8vO(A  +  O.O2)-°  4] 


(20) 


(21) 


where  (p  is  the  wind  direction  (degrees)  relative  to  the  antenna  boresight.  The  wind  speed  factor  is 


where  the  exponent  q*  =  1.1  (2  +  0.02)  -°  4 .  Then,  the  average  radar  cross  section  per  unit  area  for 
horizontal  polarization  is 

<4=  101og(3.9  x  10-6A/ 4A,AuAw)  (dB)  (23) 


where  the  subscript  hh  indicates  horizontal  polarization  on  transmit  and  receive.  (Cross-polarization 
clutter  levels  are  not  considered.)  For  vertical  polarization  and  frequencies  of  3000  MHz  and  greater 

=  1.05  ln(/ta*  + 0.02)  +  1.09  ln(2)  +  1.271n(y+  0.0001)  +  9.7  (24) 

(dB) 

where  In  is  the  natural  logarithm.  For  vertical  polarization  and  frequencies  less  than  3000  MHz 

=  (4,  -  1.73!n(A„y  +  0.02)  +  3.76  ln(2)  +  2.46  ln(V-  +  0.0001)  +  22.2  (25) 

(dB) 
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For  circular  polarization,  o°  Is  equal  to  the  larger  of  o$*  or  minus  6  dB.  The  clutter  at  the  maximum 
range  is  determined  by  linear  interpolation.  If  the  maximum  range  exceeds  the  range  for  which  grazing 
angles  can  be  calculated,  an  attenuation  rate  is  applied  to  the  last  calculated  value  of  o^at  the  limiting 
grazing  angle  Vv  The  attenuation  rate  is  adapted  from  the  standard  propagation  model:5 

a  =  {92.516  -  [8608.7593  -  (A  -  20.2663)Y/2}  (dB/km)  (26) 


where  A  is  the  scaled  evaporation  duct  height,  end  a  is  constrained  to  be  no  less  than  0.0009.  The 
scaled  evaporation  duct  height  is 

A  =  6  ZN  (m)  (27) 


which  is  constrained  to  be  23.3  meters  or  less.  The  attenuation  rate  is  also  scaled  by 

P  =  a  Rn  (dB/km)  (28) 


Zn  and  Rn  are  height  and  range  scaling  factors  found  from 

ZN  =  2.214  lO'3/2'3 
Rn  =  4.705  10'2/1'3 


The  normalized  radar  clutter  cross  section  at  the  maximum  range  is  then 


Ormu  =  <$l- 2/3  Ar 


(dB) 


(30) 


where  Ar  is  the  range  difference  between  rmax  and  the  range  of  Vv  This  model  for  normalized  radar 
clutter  cross  section  does  not  take  ocean  swell  into  account. 

The  clutter  cross  section  is 


oc  =  o°  +  Ac  (dBsm)  (31) 

where  <f  is  normalized  radar  cross  section  of  the  sea  surface  and  Ac  is  the  area  of  ocean  surface 
illuminated  by  the  radar: 


Ac  =  10  log  (1890  BWh  iC  ■ 


(dBsm) 


(32) 


where  BWh  is  horizonta,  beamwidth  in  degrees,  rc  is  compressed  pulse  length  in  microseconds,  and  r  is 
range  in  kilometers.  If  the  radar  does  not  use  pulse  compression  techniques,  then  the  normal  pulse  length 
in  microseconds  is  used  in  determining  the  area  of  the  radar  resolution  cell. 
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4.0  TEST  CASES  FOR  CALCULATION  OF  AVERAGE  SEA  CLUTTER 


This  section  lists  input  data  and  output  data  to  verify  the  proper  implementatation  of  the  computer 
program  listed  in  the  appendix.  Intermediate  results  are  also  included  to  aid  in  isolating  errors.  Table  3 
lists  the  input  data  for  the  first  test  case.  This  test  case  is  typical  of  clutter  under  standard  propagation 
conditions  for  an  X-band  radar  system.  Table  4  shows  intermediate  calculations  of  elevation  angle, 
grazing  angle  at  the  surface,  normalized  radar  sea  clutter,  and  radar  resolution  cell  versus  range.  Table  S 
shows  the  final  results  of  sea  clutter  versus  range.  If  the  program  is  properly  implemented,  the  final  results 
should  agree  within  ±1  in  the  least  significant  digit. 

Table  3.  Input  data  for  Test  Case  1. 


Input 

Description 

5.0 

1  surface  wind  speed  in  m/s 

0.0 

1  wind  dir.  rel.  to  boresight  (deg) 

0.0 

1  evap  duct  ht  in  m 

4 

t  #  of  levels  in  profile 

0.0 

350.0 

I  ht.  &  M  arrays  in  m  &  M 

10.0 

351.0 

1000.0 

468.0 

5000.0 

940.0 

30.0 

!  radar  antenna  ht  (m) 

9000.0 

1  radar  freq  (MHz) 

H 

1  antenna  polarization 

2.5 

I  hor.  beam  width  (deg) 

1.0 

1  compressed  pulse  length  (>is) 

50.0 

!  max  range  (km) 
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Table  4.  Intermediate  results  for  Test  Case  1. 


Range, 
r  (km) 

Elevation  Angle, 
a  (rad) 

Grazing  Angle, 
V'  (rad) 

Normalized  Radar 
Sea  Clutter, 
o°(dB) 

Radar  Resolution 
Cell, 

Ac  (dBsm) 

Radar  Sea 
Clutter, 
oc  (dBsm) 

0.17452 

-0.17193 

0.17191 

-39.82479 

29.16253 

-10.66226 

0.23234 

-0.12912 

0.12909 

-39.98860 

30.40524 

-9.58336 

0.34773 

-0.08630 

0.08626 

-40.36736 

32.15649 

-8.21087 

0.38600 

-0.07774 

0.07769 

-40.49141 

32.60989 

-7.88152 

0.43380 

-0.06918 

0.06912 

-40.64383 

33.11694 

-7.52689 

0.49523 

-0.06061 

0.06055 

-40.83755 

33.69205 

-7.14550 

0.57673 

-0.05205 

0.05198 

-41.09835 

34.35378 

-6.74457 

0.69049 

-0.04349 

0.04340 

-41.48647 

35.13559 

-6.35088 

0.86033 

-0.03492 

0.03482 

-42.17085 

36.09069 

-6.08016 

1.14118 

-0.02636 

0.02622 

-43.71303 

37.31754 

-6.39549 

1.69574 

-0.01780 

0.01759 

-47.98866 

39.03761 

-8.95105 

1.78255 

-0.01694 

0.01672 

-48.71342 

39.25444 

-9.45898 

1.87883 

-0.01608 

0.01585 

-49.51404 

39.48289 

-10.03115 

1.98618 

-0.01523 

0.01498 

-50.39718 

39.72420 

-10.67298 

2.10665 

-0.01437 

0.01411 

-51.37030 

39.97995 

-11.39035 

2.24286 

-0.01352 

0.01324 

-52.44196 

40.25204 

-12.18991 

2.39813 

-0.01266 

0.01236 

-53.62233 

40-54275 

-13.07958 

2.57681 

-0.01180 

0.01148 

-54.92403 

40.85485 

-14.06919 

2.78473 

-0.01095 

0.01060 

-56.36315 

41.19185 

-15.17130 

3.02988 

-0.01009 

0.00971 

-57.96095 

41.55827 

-16.40268 

3.32345 

-0.00923 

0.00882 

-59.74651 

41.95991 

-17.78660 

3.68186 

-0.00838 

0.00792 

-61.76105 

42.40469 

-19.35636 

4.13012 

-0.00752 

0.00701 

-64.06571 

42.90364 

-21.16207 

4.70864 

-0.00667 

0.00608 

-66.75675 

43.47298 

-23.28378 

5.48822 

-0.00581 

0.00512 

-69.99928 

44.13834 

-25-86095 

6.60798 

-0.00495 

0.00413 

-74.11634 

44.94470 

-29.17164 

8.39934 

-0.00410 

0.00305 

-79.90308 

45.98647 

-33.91661 

12.06655 

-0.00324 

0.00173 

-90.68350 

47.55985 

-43.12365 
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Table  5.  Final  values  for  Test  Case  1. 


Range, 
r  (km) 

Average  Radar 

Sea  Clutter, 

°c  (dBsm) 

0.17 

-10.66 

0.23 

-9.58 

0.35 

-8.21 

0.39 

-7.88 

0.43 

-7.53 

0.50 

-7.15 

0.58 

-6.74 

0.69 

-6.35 

0.86 

-6.08 

1.14 

-6.40 

1.70 

-8.95 

1.78 

-9.46 

1.88 

-10.03 

1.99 

-10.67 

2.11 

-11.39 

2.24 

-12.19 

2.40 

-13.08 

2.58 

-14.07 

2.78 

-15.17 

3.03 

-16.40 

3.32 

-17.79 

3.68 

-19.36 

4.13 

-21.16 

4.71 

-23.28 

5.^9 

-25.86 

6.61 

-29.17 

8.40 

-33.92 

12.07 

-43.12 

50.00 

-183.45 

Table  6  lists  the  input  data  for  Test  Case  2.  This  test  case  varies  from  Test  Case  1  in  that  the  relative 
wind  is  a  crosswind  and  the  propagation  conditions  are  that  of  a  13-meter  evaporation  duct.  Table  7 
shows  intermediate  calculations  of  elevation  angle,  grazing  angle  at  the  surface,  normalized  radar  sea 
clutter,  and  radar  resolution  cell  versus  range.  Table  8  shows  the  final  results  of  sea  clutter  versus  range. 
If  the  program  is  properly  implemented,  the  final  results  should  agree  within  ±1  in  the  least  significant 
digit. 
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Table  6.  Input  data  for  Test  Case  2. 


Input 

Description 

5.0 

1  surface  wind  speed  in  m/s 

90.0 

1  wind  dir.  rel.  to  boresight  (deg) 

13.0 

1  evap  duct  ht  in  m 

4 

1  #  of  levels  in  profile 

0.0 

350.0 

1  ht.  &  M  arrays  in  m  Sl  M 

10.0 

351.0 

1000.0 

468.0 

5000.0 

940.0 

30.0 

1  radar  antenna  ht  (m) 

9000.0 

1  radar  freq  (MHz) 

H 

1  antenna  polarization 

2.5 

1  hor.  beam  width  (deg) 

1.0 

1  compressed  pulse  length  (jx s) 

50.0 

1  max  range  (km) 
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Table  7.  Intermediate  results  for  Test  Case  2. 


Range, 
t  (km) 

Elevation  Angle, 
a  (rad) 

Grazing  Angle, 
xp  (rad) 

Normalized  Radar 
Sea  Clutter, 
o»(dB) 

Radar  Resolution 
Cell. 

Ae  (dBsm) 

Radar  Sea 
Clutter, 
°t  (dBsm) 

0.17330 

-0.17277 

0.17286 

-41.27011 

29.13188 

-12.13823 

0.23212 

-0.12949 

0.12961 

-41.77278 

30.40117 

-11.37162 

0.34837 

-0.08621 

0.08639 

-42.49220 

32.16438 

-10.32782 

0.38675 

-0.07755 

0.07776 

-42.68478 

32.61832 

-10.06646 

0.43536 

-0.06890 

0.06913 

-42.90594 

33.13248 

-9.77346 

0.49805 

-0.06024 

0.06050 

-43.16871 

33.71676 

-9.45195 

0.58163 

-0.05158 

0.05189 

-43.49904 

34.39051 

-9.10853 

0.69889 

-0.04293 

0.04330 

-43.95747 

35.18808 

-8.76938 

0.87557 

-0.03427 

0.03473 

-44.71241 

36.16693 

-8.54547 

1.17170 

-0.02562 

0.02623 

-46.30869 

37.43217 

-8.87651 

1.77088 

-0.01696 

0.01787 

-50.42794 

39.22591 

-11.20202 

1.86633 

-0.01610 

0.01705 

-51.09497 

39.45390 

-11.64107 

1.97272 

-0.01523 

0.01624 

-51.82214 

39.69468 

-12.12746 

2.09198 

-0.01436 

0.01543 

-52.61243 

39.94959 

-12.66284 

2.22664 

-0.01350 

0.01463 

-53.46852 

40.22051 

-13.24801 

2.37983 

-0.01263 

0.01383 

-54.39278 

40.50948 

-13.88329 

2.55577 

-0.01177 

0.01305 

-55.38715 

40.81924 

-14.56791 

2.75992 

-0.01090 

0.01227 

-56.45305 

41.15298 

-15.30007 

2.99960 

-0.01004 

0.01151 

-57.59114 

41.51465 

-16.07649 

3.28506 

-0.00917 

0.01077 

-58.80091 

41.90945 

-16.89146 

3.63102 

-0.00831 

0.01004 

-60.08001 

42.34431 

-17.73570 

4.05901 

-0.00744 

0.00934 

-61.42309 

42.82822 

-18.59487 

4.60268 

-0.00657 

0.00866 

-62.82013 

43.37412 

-19.44600 

5.31701 

-0.00571 

0.00802 

-64.25374 

44.00069 

-20.25305 

6.29973 

-0.00484 

0.00743 

-65.69566 

44.73724 

-20.95842 

7.74440 

-0.00398 

0.00690 

-67.10223 

45.63390 

-21.46833 

10.10779 

-0.00311 

0.00644 

-68.41035 

46.79058 

-21.61978 

14.89990 

-0.00225 

0.00607 

-69.53645 

48.47585 

-21.06060 

52.84929 

-0.00138 

0.00581 

-70.38317 

53.97441 

-16.40876 
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Table  8.  Final  values  for  Test  Case  2. 


Range, 
r  (km) 

Average  Radar 

Sea  Clutter, 

°c  (dBsm) 

0.17 

-12.14 

0.23 

-11.37 

0.35 

-10.33 

0.39 

-10.07 

0.44 

-9.77 

0.50 

-9.45 

0.58 

-9.11 

0.70 

-8.77 

0.88 

-8.55 

1.17 

-8.88 

1.77 

-11.20 

1.87 

-11.64 

1.97 

-12.13 

2.09 

-12.66 

2.23 

-13.25 

2.38 

-13.88 

2.56 

-14.57 

2.76 

-15.30 

3.00 

-16.08 

3.29 

-16.89 

3.63 

-17.74 

4.06 

-18.59 

4.60 

-19.45 

5.32 

-20.25 

6.30 

-20.96 

7.74 

-21.47 

10.11 

-21.62 

14.90 

-21.06 

50.00 

-16.76 

Table  9  lists  the  Input  data  for  Test  Case  3.  Test  Case  3  verifies  the  proper  operation  of  the  program 
at  a  lower  frequency  under  propagation  conditions  defined  by  an  effective  earth  radius  factor  of  1-2 
(dM/di  =  0.13  M/m).  Table  10  shows  intermediate  calculations  of  elevation  angle,  grazing  angle  at  the 
surface,  normalized  radar  sea  clutter,  and  radar  resolution  cell  versus  range.  Table  11  shows  the  final 
results  of  sea  clutter  versus  range.  If  the  program  is  properly  implemented,  the  final  results  should  agree 
within  ±  1  in  the  least  significant  digit. 
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Table  9.  Input  data  for  Test  Case  3. 


Input 

Description 

5.0 

1  surface  wind  speed  in  m/s 

0.0 

t  wind  dir.  rel.  to  boresight  (deg) 

0.0 

I  evap  duct  ht  in  m 

3 

t  #  of  levels  in  profile 

0.0 

350.0 

1  ht.  A  M  arrays  in  m  &  M 

10.0 

480.0 

1000.0 

1000.0 

5000.0 

30.0 

1  radar  antenna  ht  (m) 

1000.0 

1  radar  freq  (MHz) 

H 

1  antenna  polarization 

2.4 

1  hor.  beam  width  (deg) 

1.0 

t  compressed  pulse  length  (ns) 

50.0 

t  max  range  (km) 
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Table  10.  Intermediate  results  for  Test  Case  3. 


Range, 
r  (km) 

Elevation  Angle, 
o  (rad) 

Grazing  Ar.gle, 

\p  (rad) 

Normalized  Radar 
Sea  Clutter, 
o°(dB) 

Radar  Resolution 
Cell, 

Ac(dBsm) 

Radar  Sea 
Clutter, 
oc  (dBsm) 

0.17188 

-0.17453 

0.17453 

-47.58787 

28.91886 

-18.66901 

0.43307 

-0.06925 

0.06925 

-58.02119 

32.93227 

-25.08892 

0.69426 

-0.04317 

0.04317 

-66.57854 

34.98193 

-31.59661 

0.95545 

-0.03134 

0.03134 

-72.60278 

36.36879 

-36.23398 

1.21664 

-0.02458 

0.02458 

-77.20615 

37.41834 

-39.78781 

1.47783 

-0.02020 

0.02020 

-80.93079 

38.26297 

-42.66782 

1.73902 

-0.01714 

0.01714 

-84.06212 

38.96977 

-45.09235 

2.00021 

-0.01487 

0.01487 

-86.76709 

39.57748 

-47.18961 

2.26140 

-0.01312 

0.01312 

-89.15154 

40.11050 

-49.04104 

2.52259 

-0.01173 

0.01173 

-91.28673 

40.58520 

-50.70153 

2.78378 

-0.01060 

0.01060 

-93.22287 

41.01308 

-52.20979 

3.04497 

-0.00966 

0.00966 

-94.99673 

41.40256 

-53.59417 

3.30616 

-0.00886 

0.00886 

-96.63602 

41.75997 

-54.87605 

3.56735 

-0.00818 

0.00818 

-98.16217 

42.09019 

-56.07198 

3.82854 

-0.00759 

0.00759 

-99.59206 

42.39706 

-57.19499 

4.08973 

-0.00707 

0.00707 

-100.93926 

42.68368 

-58.25558 

4.35092 

-0.00661 

0.00661 

-102.21488 

42.95255 

-59.26234 

4.61211 

-0.00621 

0.00621 

-103.42805 

’.20573 

-60.22232 

4.87330 

-0.00584 

0.00584 

-104.58645 

43.44497 

-61.14148 

5.13449 

-0.00551 

0.00551 

-105.696S7 

43.67171 

-62.02486 

5.39568 

-0.00521 

0.00521 

-106.76396 

43.88720 

-62.87677 

5-65688 

-0.00494 

0.00494 

-107.79343 

44.09249 

-63.70094 

5.91807 

-0.00469 

0.00469 

-108.78914 

44.28853 

-64.50061 

6.17926 

-0.00446 

0.00446 

-109.75471 

44.47609 

-65.27863 

6.44045 

-0.00424 

0.00424 

-110.69342 

44.65589 

-66.03753 

6.70164 

-0.00404 

0.00404 

-111.60809 

44.82854 

-66.77954 

6.96283 

-0.00386 

0.00386 

-112.50133 

44.99459 

-67.50674 

7.22402 

-0.00369 

0.00369 

-113.37543 

45.15452 

-68.22092 

7.48521 

-0.00352 

0.00352 

-114.23251 

45.30877 

-68.92374 

7.74640 

-0.00337 

0.00337 

-115.07449 

45.45773 

-69.61676 

3  00759 

-0.00323 

0.00323 

-115.90313 

45.60175 

-70.30138 

8.26878 

-0.00309 

0.00309 

-116.72008 

45.74115 

-70.97893 

8.52997 

-0.00297 

0.00297 

-117.52682 

45.87621 

-71.65061 

8.79116 

-0.00284 

0.00284 

-118.32481 

46.00719 

-72.31762 

9.05235 

-0.00273 

0.00273 

-119.11538 

46.13435 

-72.98103 

9.31354 

-0.00262 

0.00262 

-119.89980 

46.25788 

-73  64191 

9.57473 

-0.00251 

0.00251 

-120.67928 

46.37800 

-74.30128 

9.83592 

-0.00241 

0.00241 

-121.45501 

46.49488 

-74.96013 

10.09711 

-0.00232 

0.00232 

-122.22810 

46.60870 

-75.61940 

10.35830 

-0.00223 

0.00223 

-122.99966 

46.71962 

-76.28004 

10.61950 

-0.00214 

0.00214 

-123.77080 

46.82777 

-76.94302 

10.88069 

-0.00205 

0.00205 

-124.54256 

46.93330 

-77.60926 

11.14188 

-0.00197 

0.00197 

-125.31601 

47.03632 

-78.27969 

11.40307 

-0.00189 

0.00189 

-126.09224 

47.13695 

-78.95529 

11.92545 

-0.00175 

0.00175 

-127.65739 

.7.33148 

-80.32591 
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Table  11.  Final  values  for  Test  Case  3. 


Range, 
r  (km) 

Average  Radar 

Sea  Clutter, 
oe  (dBsm) 

0.17 

-18.67 

0.43 

-25.09 

0.69 

-31.60 

0.96 

-36.23 

1.22 

-39.79 

1.48 

-42.67 

1.74 

-45.09 

2.00 

-47.19 

2.26 

-49.04 

2.52 

-50.70 

2.78 

-52.21 

3.04 

-53.59 

3.31 

-54.88 

3.57 

-56.07 

3.83 

-57.19 

4.09 

-58.26 

4.35 

-59.26 

4.61 

-60.22 

4.87 

-61.14 

5.13 

-62.02 

5.40 

-62.88 

5.66 

-63.70 

5.92 

-64.50 

6.18 

-65.28 

6.44 

-66.04 

6.70 

-66.78 

6.96 

-67.51 

7.22 

-68.22 

7.49 

-68.92 

7.75 

-69.6? 

8.01 

-70.30 

8.27 

-70.98 

8.53 

-71.65 

8.79 

-72.32 

9.05 

-72.98 

9.31 

-73.64 

9.57 

-74.30 

9.84 

-74.96 

10.10 

-75.62 

10.36 

-76.28 

10.62 

-76.94 

10.88 

-77.61 

11.14 

-78.28 

11.40 

-78.96 

11.93 

-80.33 

50.00 

-144.79 

Table  12  lists  the  input  data  for  Test  Case  4.  Test  Case  4  tests  the  program  for  a  higher  antenna 
height,  vertical  polarization,  a  wider  horizontal  beam,  and  a  longer  compressed  pulse.  Because  of  the 
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higher  antenna,  the  maximum  range  has  been  set  to  200  km.  Table  13  shows  intermediate  calculations 
of  elevation  angle,  grazing  angle  at  the  surface,  normalized  radar  sea  clutter,  and  radar  resolution  cell 
versus  range.  Table  14  shows  the  final  results  of  sea  clutter  versus  range.  If  the  program  is  properly 
implemented,  the  final  results  should  agree  within  ±1  in  the  least  significant  digit. 

Table  12.  Input  data  for  Test  Case  4. 


Input 

Description 

5.0 

1  surface  wind  speed  in  m/s 

0.0 

1  wind  dir.  rel.  to  boresight  (deg) 

0.0 

1  evap  duct  ht  in  m 

4 

1  #  of  levels  in  profile 

0.0 

350.0 

1  ht.  &  M  arrays  in  m  St  M 

10.0 

351.0 

1000.0 

468.0 

5000.0 

940.0 

1000.0 

t  radar  antenna  ht  (m) 

9000.0 

1  radar  freq  (MHz) 

V 

1  antenna  polarization 

5.0 

1  hor.  beam  width  (deg) 

2.5 

1  compressed  pulse  length  (ps) 

1  max  range  (km) 

Table  13.  Intermediate  results  for  Te;t  Case  4. 


Range,  Elevation  Angle, 
r  (km)  a  (rad) 

Grazing  Angle, 
xp  (rad) 

Normalized  Radar  Radar  Resolution  Radar  Sea 

Sea  Clutter,  Cell,  Clutter, 

(dB)  ac  (dBsm)  °c  (dBsm) 

6.01684 

-0.16655 

0.16584 

-34.76330 

50. scan 

15.79501 

7.91571 

-0.12680 

0.12586 

-35.28127 

51.74952 

16.46825 

11.57980 

-0.08704 

0.08567 

-36.13678 

53.40163 

17.26485 

12.76520 

-0.07909 

0.07758 

-36.38067 

53.82490 

17.44423 

14.22467 

-0.07114 

0.06946 

-36.66499 

54.29504 

17.63005 

16.06660 

-0.06319 

0.06129 

-37.00522 

54.82386 

17.81864 

18.46782 

-0.05524 

0.05305 

-37.42987 

55.42878 

17.i  7891 

21.73692 

-0.04729 

0.04471 

-38.00035 

56.13660 

18.13626 

26.47287 

-0.03934 

0.03620 

-38.87930 

56.99263 

18.11333 

34.04010 

-0.03139 

0.02736 

-40.62358 

58.08453 

17.46095 

48.62500 

-0.02344 

0.01768 

-45.67751 

59.63322 

13.95571 

50.93113 

-0.02264 

0.01661 

-46.65054 

59.83446 

13.18392 

53.51353 

-0.02184 

0.01551 

-47.77897 

60.04926 

12.27029 

56.43565 

-0.02105 

0.01436 

-49.09777 

60.28016 

11.18238 

59. '’8544 

-0.02025 

0.01317 

-50.65598 

60.53058 

9.87460 

63.69080 

-0.01946 

0.01191 

-52.52724 

60.80539 

8.27815 

68.34922 

-0.01866 

0.01056 

-54.83195 

61.11196 

6.28001 

74.09347 

-0.01787 

0.00909 

-57.79104 

61.46242 

3.67138 

81.56604 

-0.01707 

0.00740 

-61.88675 

61.87971 

-0.00704 

92.36490 

-0.01628 

0.00532 

-68.56870 

62.41969 

-6.14901 

115.43026 

-0.01548 

0.00173 

-91.32991 

63.38782 

-27.94209 
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Table  14.  Final  values  for  Test  Case  4. 


Range, 
r  (km) 

Average  Radar 

Sea  Clutter, 

°c  (dBsm) 

6.0 

15.8 

7.9 

16.5 

11.6 

17.3 

12.8 

17.4 

14.2 

17.6 

16.1 

17.8 

18.5 

18.0 

21.7 

18.1 

26.5 

18.1 

34.0 

17.5 

48.6 

14.0 

50.9 

13.2 

53.5 

12.3 

56.4 

11.2 

59.8 

9.9 

63.7 

8.3 

68.4 

6.3 

74.1 

3.7 

81.6 

0.0 

92.4 

-6.2 

115.4 

-27.9 

200.0 

-352.2 
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Appendix 


FORTRAN  PROGRAM  TO  CALCULATE  AVERAGE  RADAR 
CLUTTER  CROSS  SECTION  OF  THE  SEA 


ccccccccccccccccccccccccccccccccccccccceccccccccccccccccccccccccccccccc 

c 

c 

c 

c  Name : 
c 

C  PROGRAM  MAIN 

C 

c  Purpose : 
c 

c  This  program  generates  values  of  average  radar  clutter  cross 

c  section  of  the  sea  (dB)  versus  range  (km) . 

c 

c  Inputs: 
c 

c  cpl,  delta,  freq,  hbw,  height(*),  Ht,  im&x,  levels.  Munits(»), 

c  phi,  polar,  r("),  rmax,  sigmacC),  Vw 

c 

c  Outputs: 
c 

c  r<*) ,  sigmac(*) 

c 

c  Routines  called: 
c 

c  clutr,  system 

c 

c  Glossary: 
c 

c  cpl 

c  delta 

c  freq 

c  hbw 

c  height (*) 

c  Ht 

c  imax 

c  levels 

c  Munits(*) 

c  phi 

c 

c  polar 

c 


-  compressed  pulse  length  (us) 

-  evaporation  duct  height  (0  to  40  m) 

-  radar  frequency  (100  to  20000  MHz) 

-  radar  horizontal  beamwidth  (deg) 

-  array  of  height  values  for  M-profile  (m) 

-  radar  antenna  height  (m) 

-  number  of  elements  in  r  and  sigmac  arrays 

-  number  of  elements  in  neight  &  Munits  arrays 

-  array  of  M-values  (M) 

-  surface  wind  direction  relative  to  antenna 
boresight  (0  to  180  deg) 

-  antenna  polarization 

H  -  horizontal 
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r<») 

rmax 

sigmac(*) 

Vw 


V  -  vertical 
C  -  circular 

array  of  range  values  (km) 
maximum  range  (km) 
array  of  average  sea  res  (dB) 
wind  speed  (m/s) 


Method: 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c 

PROGRAM  MAIN 


Input  data  are  read  from  text  files  and  passed  to  the  clutter 
subroutine  for  calculations.  Average  radar  clutter  cross 
section  versus  range  is  returned  and  printed  out.  Prompts 
for  file  names  are  written  to  'stderr'  and  output  is  written 
to  'stdout'  to  facilitate  redirection  of  output. 


c 

c 

DIMENSION  height(30),  sigmac(50),  r(50) 
REAL  •  4  Munits (30) 
c 

CHARACTER  *  12  filnam 
CHARACTER  *  1  polar 
INTEGER  •  2  levels 


COMMON  /ENVIN/  Vw,  phi,  delta,  levels,  height,  Munits 
COMMON  /SYSIN/  Ht,  freq,  polar,  hbw,  cpl,  rmax 
c 

CALL  SYSTEM  ('Is  E*  1>&2' //char(O) ) 

WRITE  (7 ,'(//, "Enter  input  environment  file  name")') 

READ  (5. ' (al2) ' )  filnam 
OPEN  (10, FILE-f ilnam) 
c 

READ  (10,  '(f4.1)')  Vw  !  surface  wind  speed  (m/s) 

READ  (10,  '(f3,0)')  phi  !  wind  dir.  rel.  to 

bores ight 

READ  (10.  ' ( f 3 . 0) ' )  delta  !  evap  duct  ht  (m) 

READ  (10,  '(12)')  levels 
DO  i — 1 ,  levels 

READ  (10.  ' (2fl0. 1) ' )  height(i),  Munits(i)  !  ht  k  m-units 
END  DO 
CLOSE  (10) 


C 


CALL  SYSTEM  ( ' Is  S»  1>&2' //char(O) ) 

WRITE  (7 ,'(//, "Enter  input  system  file  name")') 
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READ  (5 , ' (al2) ' )  filnam 
OPEN  (lO.FILE-filnam) 


c 


c 

c 

c 


c 


READ  (10, 
READ  (10. 
READ  (10, 
READ  (10, 
READ  (10, 
READ  (10. 
CLOSE(IO) 


'(flO.l)')  Ht 
'(flO.l)')  freq 
'(al)')  polar 
'(flO.l)')  hbw 
' (flO.l) ' )  cpl 
' (£4.0) ' )  rmax 


CALL  clutr  (sigmac,  r,  imax) 


!  radar  antenna  ht  (m) 

!  radar  freq  (MHz) 

!  antenna  polarization 
!  hor.  beam  width  (deg) 

!  compr.  pulse  length  (us) 
!  max  range  (km) 


DO  i  -  1,  imax 

WRITE  (6,  ' (2f 10. 2) ' )  r(i),  sigmac(i) 
END  DO 

STOP 

END 


ccccccccccccccccccccccccccccccocccccccccccccccccccccccccccccccccccccccc 

c 

c 

c 

c  Name : 
c 

c  SUBROUTINE  clutr  (sigmac,  r,  imax) 

c 

c  Purpose: 
c 

c  Calculates  average  radar  cross  section  of  the  sea  versus 

c  range . 

c 

c  Inputs: 
c 

c  alpha(*),  cpl,  delta,  freq,  hbw,  imax,  psi(*),  r(*), 

c  rmax,  sigma0(*) 

c 

c  Outputs: 
c 

c  sigmac(»),  r(*) 

c 

c  Calling  routines: 
c 

c  MAIN 

c 

c  Routines  called: 
c 

c  AL0G10,  FLOAT,  getk,  IFIX,  mprof,  nrcs,  psir,  psirae,  SORT 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Glossary: 

alph&( • ) 

cpl 

delta 

freq 

hbw 

imax 

P«i(*) 

r{») 

rmax 

sigmac(*) 

sigmaOC) 


-  depression  angle  array,  deg. 

-  radar  receiver  compressed  pulse  length,  microsec 

-  evaporation  duct  height,  0-40  meters 

-  radar  system  frequency,  100  -  20000  MHz 

-  radar  antenna  horizontal  beam  width,  deg 

-  number  of  elements  in  sigmac  array 

-  grazing  angle  array 

-  range  array  corresponding  to  sigmac  (km) 

-  maximum  range,  km 

-  average  radar  cross  section  of  the  sea  dB 

-  normalized  radar  reflectivity  of  the  sea  surface, 
sq.  meters  per  sq.  meter 


Method: 


A  grazing  angle  model  is  selected,  depending  upon 
frequency,  to  determine  grazing  angle  versus  range.  These 
grazing  angles  are  input  to  the  routine  that  calculates 
normalized  sea  clutter  cross  section.  Normalized  clutter 
is  converted  to  average  clutter  by  multiplying  by  the 
radar's  clutter  resolution  cell  at  each  range.  Beyond 
the  range  of  the  minimus1  grazing  angle,  an  attenuation 
rate  dependent  upon  evaporation  duct  conditions  is  applied 
to  the  sea  surface  reflectivity  calculated  at  the  range  of 
the  minimum  grazing  angle. 


c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c 

SUBROUTINE  clutr  (sigmac,  r,  imax) 
c 

DIMENSION  height(30) 

REAL  •  4  Munits(30) 

DIMENSION  alpha (SO) ,  psi(50),  r(50) ,  sigma0(50),  sigmac(50) 
DIMENSION  dMdh( 32) ,  hmrs(32) 
c 

CHARACTER  *  1  polar 
INTEGER  *  2  levels 
c 

COMMON  /ENVIN/  Vw,  Dhi ,  delta,  levels,  height,  Munits 
COMMON  /SYSIN/  Ht,  freq,  polar,  hbw,  cpl,  rmax 
c 

c  If  freq  >■  2  GHz,  then  evaporation  duct  effects  dominate 

c  determination  of  grazing  angle.  If  freq  <  2  GHz,  then 
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c  effective  earth  radius  effects  are  assumed  to  dominate  the 

c  determination  of  grazing  angle. 

c 

IF  (freq  .GE.  2000.)  THEN 

CALL  psir(delta,  Ht,  alpha,  psi,  r,  imax) 

ELSE 

CALL  mprof  (height,  Munits,  Ht,  levels,  ALPHAC,  DMDH,  HMRS, 

2  SBDHT,  NTOT) 

CALL  getk  (alphac,  dMdh,  hmrs,  ntot,  Ht,  RK) 

CALL  psirae  (rk,  Ht.  rmax,  alpha,  psi,  r,  imax) 

END  IF 
c 

c  invert  r,  psi,  and  alpha  arrays  so  range  increases  with 

index 

c 

limit-  I FIX ( FLOAT ( imax) /2 . ) 

DO  i-1,  limit 
rng  -  r ( i ) 
angle  -  alpha(i) 
angl  -  psi(i) 
r(i)  -  r(imax+l-i) 
psi(i)  -  psi (imax+l-i) 
alpha (i)  -  alpha ( imax+l-i ) 
r (imax+l-i)  -  rng 
psi (imax+l-i)  -  angl 
alpha (imax+l-i)  -  angle 
END  DO 
c 

c  determine  reflectivity  from  grazing  angle 

c 

DO  i  -  1,  imax 

CALL  nrcs  (Vw,  phi,  freq,  polar,  psi(i),  sigmaO(i)) 

END  DO 
C 

c  calculate  average  clutter  cross  section  from  normalized 

c  reflectivity  and  clutter  resolution  cell  size 

c 

acs  -  10.*AL0C10(1890.»hbw»cpl) 

DO  i  -  1 ,  imax 

ac  -  acs  +  10.*AL0G10(r(i)) 
sigmac(i)  -  sigmaO(i)  +  ac 

write  (0, ' (6fl2 . 5) ' )  r(i),  alpha(i),  psi(i),  sigmaO(i),  ac, 

2  sigmac(i) 

END  DO 
C 

c  beyond  minimum  grazing  angle,  apply  an  attenuation  rate 

c  to  the  reflectivity 

c 

IF  (r(imax)  .LT,  rmax)  THEN 
reflim  •  sigmaO(imax) 
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rlim  -  r(imax) 

zfac  -  (freq/9600. ) •* (2 . /3 . ) 
rfac  -  (freq/9600. )**(l./3. ) 
del  -  delta  *  zfac 
IF  (del  .GT.  23.3)  del  -  23.3 

atten  -  92.516  -  SQRT(8608. 7S93-(del-20. 2663) *»2) 
ir  (atten  .LT.  0.0009)  atten  -  0.0009 
atten  -  atten  *  rfac 

sigma  -  reflia  -  2.  *  atten  *  (rmax  -  rlim) 
ac  -  acs  +  10. *AL0010(raax) 
iaax  “  laax  +  1 
r(ioax)  -  raax 
sigmac(imax)  -  sigma  +  ac 
ELSE  IF  (r(imax)  .GT.  raax)  THEN 
c 

c  use  linear  interpolation  to  find  sigma_c  at  rmax 

c 

1-1 

DO  WHILE  (rmax  .GT.  r(i)) 
i  -  i  ♦  1 
END  DO 

dr  -  (roax-r ( i-1) )  /  (r(i)-r(i-l)) 
da  -  sigmac(i)  -  sigmac(i-l) 
temp  -  ds  •  dr  ♦  sigaac(i-l) 
i  max  -  i 
r(imax)  «  rmax 
sigmac(imax)  -  temp 
END  IF 
c 

RETURN 

END 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c  Name : 
c 

c  SUBROUTINE  downgo  (alpha,  ixmtr,  r) 

c 

c  Purpose: 

c 

c  Ray  trace  for  downgoing  rays,  alpha  <  0 

c 

c  Inputs: 

c 


c 

c 


alpha,  rhl(*),  rml(*),  ixmtr 


A-6 


c  Outputs: 
c 

c  r 

c 

c  Calling  routines: 
c 

c  psir 

c 

c  Routines  called: 
c 

c  SORT 

c 

c  Glossary: 
c 

c  alpha  -  launch  angle  in  radians 

c  ixotr  -  index  of  level  in  profile  for  transmitter 

c  rhl(*)  -  height  (m) 

c  rml(*)  -  modified  refractivity  (M) 

c  r  range  (km) 

c 

c  Method: 
c 

c  A  downgoing  ray  is  traced  through  the  successive  layers  of 

c  the  M  profile  to  determine  the  range  at  which  it  strikes  the 

c  surface, 

c 
c 

cccccccoccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

SUBROUTINE  downgo  (alpha,  ixmtr,  r) 
c 

DIMENSION  rhl(10),  rml(10) 
c 

COMMON  /RAYTRACE/  rhl,  rml 
c 

i  ■  ixmtr 

DO  WHILE  (i  .GT.  1) 

dh  -  rhl(i)  -  rhl(i-l) 
dm  ■  rml(i)  -  rml(i-l) 
dmdh  -  dm/dh 

IF  (dmdh  .EQ.  0.)  dmdh  -  le-8 

alphap  ■  -SQRT(alpha*»2  +  0.000002*dmdh*(rhl(i-l)-rhl(i) )) 
rp  -  r  +  (alphap-alpha)/(dmdh*0.001) 
alpha  -  alphap 
r  -  rp 
i  -  i  -  1 
END  DO 
C 

RETURN 
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END 


C 

c  Subroutine  ducts 
c 

c  DUCTS  builds  an  array  containing  the  top,  bottom,  and 
c  minimum  refractivity  of  all  the  major  ducts  in  the 
c  atnosphere  refractivity  profile, 
c 

c  Variable: 
c  dct 

c 
c 
c 

c  lvls 

c  ndcts 

c 

c  nq 

c  rht 

c  rmu 

c 
c 
c 

SUBROUTINE  ducts (rmu, rht , lvls , DCT, NDCTS) 

c 

real *4  dct .delu.delh, deltu.hbot , htop, rht (32) , reu(32) 
integer*2  lvls. ibot, lduct, lend, iq, itop, ndcts, nq 
dimension  dct(3,8) 
c 

c  Locate  all  major  ducts 

nq-0 

iq-3*ndcts 
Itop-lvla 
lend -ndcts 
ndcts-0 

DO  iduct-l.iend 
c 

c  Look  for  top  of  next  duct 

1010  continue 

htop-rht(itop) 
if (itop.eq.l)  go  to  1060 
ibot-i top-1 

if (rmu(itop) . le.rmu(ibot) )  go  to  1020 
itop-itop-1 
go  to  1010 
c 

c  Look  for  bottom  of  the  duct 

1020  continue 

hbot-rht(ibot) 

if (rmu(ibot) . It .rmu(itop) )  go  to  1030 
if (ibot.eq. 1)  go  to  1040 


Description: 

3,*  duct  parameters  array. 

1, n  bottom  of  duct  'n' ,  meters. 

2, n  top  of  duct  'n' ,  meters. 

3, n  minimum  refractivity  of  duct  'n',  M-units. 
Number  of  refractivity  level  in  rmu,  rhta. 

in:  the  maximum  number  of  ducts  allowed, 
out:  the  number  of  ducts  found. 

Duct  counter. 

Height  array,  meters. 

Modified  refractivity,  M-unit  array,  elements 
correspond  to  like-number  elements  of  rht  array. 
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ibot-ibot-l 
go  to  1020 
c 

c  Calculate  bottom  of  duct  using  linear  interpolation 

1030  continue 

delu-rmu( ibot+1 ) -rmu ( ibot ) 
delh-rh  t ( ibot+1 ) -rht ( ibot ) 
de 1 tu-rmu ( i top ) -rmu ( i bo t ) 
if (delu.lt. 0.01)  go  to  1040 
hbot-rht(ibot)  +  deltu*delh/delu 
c 

c  Store  duct  parameters  in  array  dct 

1040  continue 

amu-rmu(itop) 
call  push (dct, iq.nq.amu) 
call  push(dct, iq.nq.htop) 
call  push(dct, iq.nq.hbot) 
ndcts-iduct 
i top- ibot 
END  DO 
c 

1060  continue 
RETURN 
END 
C 

c  Subroutine  getk 
c 

c  Subroutine  GETK  is  used  to  determine  the  effective  earth  radius 
c  factor  k.  Getk  accomplishes  this  by  tracing  a  ray  from  the  trans- 
c  mi  tier  height  to  200  NMi  (370  km).  The  ray  launch  angle  is  0  deg. 
c  if  no  surface-based  duct  exists,  or  alphac,  the  critical  angle  if 
c  one  does . 
c 

c  Variable: 
c  alphac 
c 

c  aO 

c  al 

c  deld 

c  delh 

c  delM 

c  delmdh 

c  dMdh 

c  hlast 

c  hers 

c  ntot 

c  rdeld 

c  rmax 

c  rng 

c  rk 


Description: 

Critical  angle  necessary  to  escape  duct.  If  alphac 
•  0  then  no  surface-based  duct  exists. 

Initial  ray  launch  angle,  radians. 

Ray  angle  at  top  of  layer,  radians. 

Range  difference,  km. 

Height  difference,  meters. 

M-unit  difference. 

M-unit  gradient. 

M-unit  gradient  array. 

Height  at  370  km. 

Array  of  height  elements,  in  meters. 

Maximum  number  of  elements  in  hmrs  and  dMdh  arrays. 
Range  incremented  in  ray  trace. 

Maximum  range  for  ray  trace  -  370  km. 

Range,  km. 

Effective  earth  radius  factor. 
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c  xatr  Transmitter  height  in  meters, 

c 

SUBROUTINE  getk(alphac,  dMdh,  hors,  ntot,  xmtr,  RK) 
c 

real*4  alphac,  aO,  al.  deld,  delh,  delM.  delmdh,  dMdh (32) 
real*4  blast,  hmrs(32),  rdeld,  rmax,  rng,  rk,  xmtr 
integer*2  ntot,  i 
c 

rmax  -  370.0 
h  -  xmtr 
rng  -  0.0 
aO  -  alphac 

c  Loop  to  trace  ray  through  the  atmospheric  layers. 

DO  i-2,ntot-l 

dels  -  (hmrs(i+l)  -  h)»ciMdh(i)«1.0E-3 
al  -  SQRT(aO*aO  +  2.0«dela) 
deld  -  (al  -aO)/dMdh(l) 
rdeld  •  rng  +  deld 
IF(rdeld  .OT.  rmax)  GOTO  1000 
aO  -  al 
h  -  hmrs(i-fl) 
rng  »  rdeld 
END  DO 
i  -  ntot 
1000  continue 

c  Ray  trace  in  final  layer  to  range  rmax. 

deld  -  rmax  -  rng 
al  -  aO  ♦  dMdh(i)  •  deld 
delkl  -  (al*al  -  a0»a0)»0.5 
delh  -  1000 . 0*delM/dMdh ( i ) 
hlast  -  hmrs(i)  ♦  delh 

c  Determine  the  equivalent  single-gradient  atmosphere  that 

c  would  be  required  to  trace  a  ray  launched  at  alphac  that 

c  would  arrive  at  height  -  hlast  at  a  range  of  370  km. 

delmdh  -  (-alphac) *2 .0/ rmax  +  2.0E-3*(hlast  -  xmtr) / (rmax»rmax) 
rk  -  1.0/(6371.0  •  delmdh) 

IF(rk  .CT.  5.0)  rk  -  5.0 
!F(rk  .LE.  0.5)  rk  -  0.50 
RETURN 
END 
c 

c  Subroutine  insrt 

c 

c 

c  INSRT  inserts  (or  appends)  a  new  level  into  the  M-unit  profile.  It 
c  does  this  by  locating  the  new  height  relative  to  the  existing  pro- 
c  file  heights.  If  the  new  height  is  greater  than  the  top  level,  then 
c  append  a  new  level  for  the  new  height.  If  the  new  height  is  between 
c  two  levels,  then  insert  a  new  level  for  the  new  height.  If  the  new 
c  height  is  equal  to  an  existing  level's  height,  do  not  add  a  new 
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o  level  for  the  new  height, 
c 

c  Variable: 
c  aau 

c  hors 

c 

c  iq 

c  hgt 

c  ipnt 

c 
c 

SUBROUTINE  inart (aau, hars, iq.hgt , ipnt) 
c 

realM  aau(32) , hars (32)  ,hgt 
integer »2  iq.ipnt 
c 

DO  i-l.iq 
ilevel-i 

IF (ABS(hgt-hars(i level) ) .LE.0.01)  go  to  1020 
IF(hars(ilevel) .CT.hgt)  go  to  1C30 
END  DO 
C 

c  Hgt  >  aau(iq) 

iq-iq+1 
ipnt-iq 

grdnt-0. 1181102 

aau (ipnt) -aau (iq-1)  +  (hgt-hmra ( iq-1) ) *grdnt 
hars (ipnt) -hgt 
go  to  1050 

c 

c  Hgt  -  hars (1 level) 

1020  continue 

ipnt-ilevel 

c  aau ( ipnt) -aau (i level) 

hars ( ipnt) -hgt 
go  to  1050 
c 

c  Hara(ilevel)  >  hgt  >  hars(ilevel-l) 

1030  continue 

c  Shift  levels  above  new  height  up  one 

DO  l-ilevel,iq 

J-iq  -  (i-ilevel) 
hars (J+l) -hors ( J ) 
aau(j+l)-aau( j) 

END  DO 
iq-iq+1 
ipnt-i level 

grdnt-(anu (ipnt+l) -aau ( ipnt-l) )/ (hars (ipnt+1) -hars (ipnt-1) ) 
aau (ipnt) -aau (ipnt-l)  +  (hgt-hars(ipnt-l) )»grdnt 
hars (ipnt) -hgt 


Description : 

Modified  refractivity  array,  M-units. 

Height  array,  waters,  each  eleaent  corresponding  to 
the  like-nuaber  aau  array  eleaent. 

Nuaber  of  levels  in  aau  and  hairs. 

Height  of  new  level  to  be  added,  meters. 

Index  pointer  to  new  level. 
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C  go  to  1050 
c 

1050  continue 
RETURN 
END 

c 

c  Subroutine  mprof 
c 

c  MPROr  nodlfies  the  M-unit  and  height  arrays  by  inserting  a  level  at 
c  the  antenna  height  using  straight  line  interpolation  (or  a  standard 
c  atmosphere  gradient)  to  calculate  its  M-unit  value.  The  new  profile 
c  is  then  used  to  locate  uny  ducts  that  might  be  contained  in  the  pro- 
c  file.  If  the  bottom  of  the  duct  is  below  the  EM  system  antenna 
c  height,  and  the  top  above  the  antenna  height,  then  a  critical  angle 
c  is  calculated  for  the  EM  system  in  the  surface-based  duct.  (It  is 
c  assumed  that  low-elevated  ducts  are  surface  ducts  if  the  EM  system 
is 

c  in  the  duct. ) 
c 

c  Variable: 
c  alphac 
duct 

c  a  mu 

c  antena 

c  antmu 

c  dcts 

c 
c 
c 

c  dMdh 

c  hbot 

c  htop 

c  height 

c  hmrs 

c 

c  lvlant 

c  lvltop 

c  Munits 

c 

c  ndcta 

c  nmax 

c  ntot 

c  rma 

c  sbdht 

c 

c  Variables  not  listed  are  temporary  variables, 

c 
c 

SUBROUTINE  mprof (height , Munits, an tens , nmax , ALPHAC , DMDH, HMRS, 

1  SBDHT, NTOT) 


Description: 

The  critical  penetration  angle  necessary  to  escape 

An  array  of  M-unit  values 
EM  system  antenna  height 

M-unit  value  at  the  EM  system  antenna  height 
24  duct  parameter  array 

1, n  bottom  of  duct  'n',  meters 

2 , n  top  of  duct  'n' ,  meters 

3, n  minimum  refractivity  of  duct  'n' ,  m-units 
M-unit  gradient  array 

Height  of  the  bottom  of  a  duct 
Height  of  the  top  of  a  duct 

Height  array  with  the  original  profile  heights 
Height  array  with  elements  corresponding  to  the  dMdh 
array  elements 
EM  system  antenna  level 

Maximum  number  of  layers  in  the  hmrs  array 
M-unit  array  with  elements  corresponding  to  the  height 
array  elements 

The  number  of  ducts  stored  in  'dcts' 

The  number  of  elements  in  the  height  and  Munit  arrays 
The  number  of  elements  in  the  dMdh  and  hmrs  arrays 
M-unit  value  at  the  minimum  on  the  duct  profile 
The  height  of  the  surface-based  duct 
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c 

c 


real*4  alphac,aou(32) , antena, dmdh( 32) ,hmrs(32) .height (30) 
real*4  Munits(30) .sbdht 
real *4  antmu.dcts.hb.ht.rma 
integer*2  lvlant, lvltop, nmax, ntot 
integer*2  ndcts 
dimension  dcts(3,8) 
c 

lvltop  »  nmax 
alphac  -  0.0 
sbdht  -  0.0 
c 

c  Copy  height  and  m-unit  arrays, 

c 

lvltop  »  nmax 
DO  i  •  1,  nmax 

hmrs(i)-height(i) 

amu(i)-Munits(i) 

END  DO 
C 

c  Insert  new  level  at  the  antenna  height, 

c 

call  lnsrt(amu.hmrs, lvltop, antena, lvlant) 
antmu-amu ( lvlant ) 
c 

c  Locate  all  major  ducts. 

ndcts-8 

call  ducts (amu, hairs,  lvltop, dots, ndcts) 
c 

c  Define  trapping  duct  parameters. 

IF(ndcts  .NE.  OJTHEH 
DO  iduct-1 .ndcts 
hb-dcts(l, iduct) 
ht«dcts(2, iduct) 
rma-dcts(3, iduct) 

IF( (antena  .CT.  hb)  .AND.  (antena  .LT.  ht))  go  to  1040 
IF(hb. It. 0.01)  go  to  1040 
END  DO 
END  IF 
c 

c  Antenna  not  inside  a  major  duct. 

go  to  1050 
c 

c  The  antenna  is  inside  a  low-level  elevated  duct 

c  or  inside  a  surface-based  duct. 

1040  continue 
sbdht  -  ht 

alphac«1.0e-3*sqrt(2.0,(antmu-rma) )  1.0e-5 

1050  continue 
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c 

o 


Delete  ell  levels  between  the  surface  and  the  antenna  level. 
DO  i  *  lvlant, lvltop 
j«i-( lvlant-2) 
hars( j)-hars(i) 
aau  ( J )  -a«u  ( 1 ) 

END  DO 
lvltop-j 
lvlant-2 

c  Calculate  the  M-unit  gradient  array, 

iend-lvltop-l 
DO  i  -  1,  lend 

delu«amu(i+l)-aau(i) 
delh-hmrs ( i+1) -hmrs { i ) 
dmdh ( i ) -1 . 0e-3*delu/delh 
END  DO 

dmdh( lvltop) -0.1181102e-3 
c 

ntot  -  lvltop 

return 

END 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c  Name: 
c 

c  SUBROUTINE  nrcs  (vw.  phi,  freq,  polar,  psi,  sigmaO) 

c 

c  Purpose : 
c 

c  Calculates  the  average  radar  clutter  cross  section  per  unit 

c  area  of  the  sea  in  sq.  o.  per  sq.  m. 

c 

c  Inputs: 
c 

c  freq,  phi,  polar,  psi,  vw 

c 

c  Outputs: 
c 

c  sigmaO 

c 

c  Calling  routines: 
c 

c  clutr 

c 

c  Routines  called: 
c 

C  ALOG,  AL0G10,  AMAX1 ,  COS,  EXP 
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Glossary: 

freq 

phi 

polar 


psi 

sigmaO 

vw 


-  frequency  (MHz) 

-  angle  between  antenna  boreslght  and  upwind 
(0  to  180  deg) 

-  antenna  polarization 

H  -  horizontal 
V  -  vertical 
C  -  circular 

-  grazing  angle  (radians) 

-  normalized  radar  cross  section  of  the  sea  (dB) 

-  wind  speed  (o/s) 


Method: 


o 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 


The  average  radar  cross  section  per  unit  area  (reflectivity) 
of  the  sea  is  calculated  from  a  deterministic  parametric  model 
developed  by  Georgia  Institute  of  Technology:  "Radar  sea 
clutter  model,"  Proc.  IEEE  International  Conference  on 
Antennas  and  Propagation,  London,  Hov.  78,  pp  6-10,  by  Horst, 
Dyer,  &  Tuley.  The  correction  for  circular  polarization 
is  made  in  accordance  with  Nathanason,  Radar  Design  Principles, 
p.  239.  Wind  and  sea  direction  and  wind  speed  and  wave 
height  are  assumed  to  be  highly  correlated:  if  this  is  not 
a  valid  assumption,  this  code  can  be  modified  to  accept 
wind  speed  and  average  wave  height  as  separate  entrys. 


c 

c 

c 

c 


c 

c 

c 


c 

c 


SUBROUTINE  nrcs  (vw,  phi,  freq,  polar,  psi,  sigmaO) 
CHARACTER*!  polar 


constants  for  clutter 

rlamda  -  300.  /  freq 
hav  -  (vw  /  8 . 87 ) **2 . 5 
cosphi  -  COS(phi/57.3) 


radar  wavelength 
average  wave  height 
cos  of  wind  direction 


interference  factor,  ai 


sigma _phi  -  (14.4  •  rlamda 

sp4-sigma_phi**4 

ai  -  sp4  /  (l.  +  sp4) 


+5.5)  *  psi  *  hav  /  rlamda 


upwind/downwind  factor,  au 
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rl4  -  (r land a  +  0.02)*»(-0.4) 
au  -  EXP(0 . 2  *  cosphi  •  (1.  -  2.8*psi)  •  rl4) 
c 

c  wind  speed  factor,  aw 

c 

qw  -  1.1  *  rl4 

aw  -  (1.0425  *  vw  /  (1.  +  vw  /  15.))**qw 
c 

c  average  radar  cross  section  per  unit  area 

c  for  horizontal  polarization,  sigma_hh 

c 

cltr  -  0.0000039  *  rlamda  •  psi**0.4  •  ai  •  au  «  aw 
sigma_hh  ■  10.  •  ALOGlO(cltr) 
sigmaO  -  sigma_hh 
c 

c  correction  factor  for  vertical  polarization 

c 

IF  ((polar  .EQ.  "V")  .OR.  (polar  . EQ.  "C"))  THEN 
IF  (freq  .GE.  3000.)  THEN 

sigma_vv  -  sigma_hh  -  1.05,AL0G(hav+0.02)  + 

1. 09 *AL0C( rlaoda) 

sigma_vv  sigma_vv  +  1.27*ALOG(psi+O.OOOi)  +  9.7 
ELSE 

sigma_vv  -  sigma_hh  -  1 . 73*AL0G(hav+0.02)  + 

3.76*AL0G(rlamda) 

sigma_vv  ■  sigma^vv  +  2. 46*AL0C(psi +0.0001)  +22.2 
END  IF 
C 

c  correction  factor  for  circular  polarization 

c 

IF  (polar  .EQ.  "C")  THEN 

sigmaO  -  AMAXl(sigma_hh, sigma_vv)  -  6 
ELSE 

sigmaO  -  sigma_vv 
END  IF 
END  IF 
c 
c 

RETURN 

END 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c 

c  Name : 
c 

c  SUBROUTINE  pair  (delta,  hxratr,  alpha,  psi,  r,  Imax) 

c 

c  Purpose; 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


This  subroutine  generates  arrays  of  ray  launch  angle,  alpha, 
grazing  angle,  psi,  and  range,  r,  to  be  used  as  inputs  to  the 
clutter  level  calculations. 

Inputs : 

delta,  hxmtr 

Outputs : 

alpha(*) ,  angle,  idelta,  imax,  ixmtr,  nmaxl,  psi(»),  r ( » ) . 
rhl(»),  rml(*) 

Calling  routines: 
clutr 

Routines  called: 

ALOG,  downgo,  EXP,  FLOAT,  SORT,  upgo 

Glossary: 

alpha(*)  -  ray  launch  angle  at  the  transmitter  (radians) 

angle  -  ray  launch  angle  at  the  transmitter  (radians) 

delta  -  evaporation  duct  height  (m) 
hxmtr  -  height  of  the  transmitter  (m) 
idelta  -  index  of  level  in  profile  for  delta 

ixmtr  -  index  of  level  in  profile  for  hxmtr 

imax  -  number  of  elements  in  alpha,  psi  and  r  arrays 

nmaxl  -  number  of  elements  in  generated  height  and  M  arrays 

(includes  radar  height  and  evap  duct  profile) 
psi(»)  -  grazing  angle  at  the  sea  surface  (radians) 
rhl(*)  -  augmented  height  array  <m) 
rml(*)  -  augmented  modified  refractivity  array  (M) 

r  -  range  (km) 

Method : 

This  subroutine  takes  the  evaporation  duct  height 
and  radar  antenna  height,  and,  using  ray  trace  techniques, 
returns  grazing  angle,  psi  (radians),  as  a  function  of  range, 
r  (m) .  Psi  is  determined  directly,  by  Snell's  law,  given  the 
radar  antenna  height,  depression  (or  elevation)  angle  at  the 
radar,  the  M-value  at  that  height,  and  the  M-value  at  the 
surface,  The  range  is  determined  by  using  small  angle 
approximations,  a  linearly  segmented  M-profile,  and  analytic 
expressions  for  the  ray  path  in  each  segment.  The  evaporation 
duct  profile  is  assumed  to  be  for  neutral  conditions  and  uses 
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c  the  Jeske  formulation.  The  evaporation  duct  profile  is 

c  extended  to  the  radar  antenna  height,  if  necessary,  by 

c  extrapolating  with  a  gradient  of  118  M/km. 

c 
c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c 

SUBROUTINE  pair  (delta,  hxatr,  alpha,  psi,  r,  imax) 
c 

DIMENSION  rhl(lO),  rml(10) 

DIMENSION  &lpha(50),  psi (50),  r(50) 

C 

COMMON  /RAYTRACE/  rhl ,  rml 


zO  0.00015  !  surface  roughness  parameter  (m) 

zl  -  6.  !  reference  height  (m) 

k  «  0.4  !  von  Karmen's  constant 

c 

c  calculate  M-profile  for  input  delta,  neutral  stability 

c 

rosfc  -  340.0 
rhl ( 1 )  -  0. 
rml(l)  -  rmsfc 

DO  i  -  2,  8  !  logarithmic 

height 

rlz  -  FLOAT (i -4) 

z  -  EXP(rlz)  !  z- 

0.135,0.368,1.0,2.7,7.4,20.1.54.6 
rhl(i)  «  z 

rml(i)  •  rmsfc  +  0.125»z  -  0.125  •  Delta  •  ALOC(z/zO) 

END  DO 
nmaxl  -  8 


c  insert  delta  into  M-profile 

c 

i  -  1 

DO  WHILE  (rhl(i)  .LT.  Delta) 
i  •  i  t  1 
END  DO 
nl  *■  i 

IF  (delta  .NE.  rhl(i))  THEN 
nmaxl  -  nmaxl  +  1 
DO  i  -  nmaxl, nl+1,-1 
rhl(I)  -  rhl (I  -  1) 
rml(l)  -  rrnld  -  1) 

END  DO 

rhl(nl)  -  delta 

rml(nl)  -  rmsfc  +  O.lES^delta  -  0. l25,delta*ALOC(delta/zO) 
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END  IF 
idelta  ■  nl 
rmmin  -  rnl(nl) 
c 

c  insert  xotr  level  into  HI  fc  Ml  arrays 

c 

IF  (hxmtr  .GT.  rhl (nmaxl))  THEN 
nmaxl  -  nmaxl  +  1 
nl  -  nmaxl 
rhl (nmaxl)  -  hxmtr 

rml(nmaxl)  -  0 . 118* (rhl (nmaxl) -rhl (nmaxl-1) )  +  rml (nmaxl-1) 
rmxmtr  -  rml (nmaxl) 

ELSEIF  (hxmtr  .  EQ.  rhl(noaxl))  THEN 
rmxmtr  «  rml (nmaxl) 
nl  -  nmaxl 
ELSE 

i  -  2 

DO  WHILE  (rhl(i)  .LT.  hxmtr) 
i  -  i  +  1 
END  DO 
nl  -  i 

IF  (hxmtr  . EQ.  rhl(i))  THEN 
rmxmtr  »  rml(i) 

ELSE 

grad  -  (rml(i)-rml(i-l) )  /  (rhl(i)  -  rhl(i-l)) 
rmxmtr  -  rml(i-l)  +  (hxmtr-rhl ( i— 1 ) ) 'grad 
nmaxl  -  nmaxl  +  1 
DO  i  -  nmaxl,  nl+1,  -1 
rhl(i)  -  rhl(i-l) 
rml(i)  -  rml(i-l) 

END  DO 

rhl(nl)  -  hxmtr 
rml(nl)  -  rmxmtr 
END  IF 
END  IF 
ixmtr  ■  nl 

IF  (hxmtr  .LT.  delta)  idelta  -  idelta  +  1 
c 

c  Determine  psi  for  -10  <  alpha  <  alphmx  deg. 

c  Alpha  is  the  launch  angle  at  the  radar;  alphmx  is 

c  the  angle  that  yields  the  longest  range  to  the  surface; 

c  Alphmn  is  max  depression  angle  -  giving  a  max  grazing 

c  angle  at  the  surface  of  approx.  10  deg. 

c 

alphmn  -  -0.1745 

radcnd  -  0.000002*(rmxmtr-rmmin) 

IF  (radcnd  . LE .  0.)  THEN 
alphmx  -  -0.000001 
ELSE 

alphmx  -  -SQRT(radcnd)  -  0.000001 
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END  ir 

psiain  -  SORT (alphmx* *2  -  0.000002* (rmxmtr-rmsfc) ) 

IP  (psiain  .It.  0.001745)  THEN 

alphmx  -  -SQRT(0. 000003+  0.000002*(rmxatr-rmsfc) ) 

END  IF 

IF  (hxmtr  .LT.  delta)  alphmx  -  -alphmx 
dalpha  -  (alphmx-alphmn)  /  200. 
i  -  0 

angle  -  alphmx 

c 

c  Use  a  variable  angle  decrement  with  smaller  steps  at  the 

c  larger  ray  launch  angles, 

c 

DO  WHILE  (angle  .GT,  alphmn) 
i  -  i  +  1 
alpha(i)  -  angle 

psi(i)  -  SQRT(alpha(i)**2  -  0.000002* (rmxmtr  -  rmsfc)) 

IF  (psi(i)  .LT.  0.0175)  THEN 
angle  *  angle  -  dalpha 
ELSEIF  (psi(i)  .LT.  0.08)  THEN 
angle  -  angle  -  dalpha* 10. 

ELSE 

angle  -  angle  -  dalpha* 50. 

END  IF 

ir  (angle  .EQ.  0.)  angle  -  0.000001 
END  DO 

imax  -  i 

c 

c  initiate  ray  trace  from  hxmtr  to  sfc 

c 

DO  i  -  1,  imax 

angle  -  alpha(i) 
rng  -  0. 

IF  (angle  .GT.  0.)  THEN 

CALL  upgo(angle,  ixatr,  idelta,  rng) 
angle  -  -alphmx 
rng  -  rng  •  2. 

END  ir 

CALL  downgo(angle,  ixmtr,  rng) 
r(i)  -  rng 
END  DO 
c 
c 

RETURN 

END 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c 
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c  Name: 
c 

c  SUBROUTING  psirae  (rk,  hxmtr,  rux,  alpha,  pal,  r,  imax) 

c 

c  Purpose : 
c 

c  This  subroutine  generates  arrays  of  ray  launch  angle,  alpha, 

c  grazing  angle,  psi,  and  range,  r,  to  be  used  as  inputs  to  the 

c  clutter  level  calculations, 

c 

c  Inputs: 
c 

c  hxmtr,  rk,  rmax 

c 

c  Outputs: 
c 

c  alpha,  psi,  r,  imax 

c 

c  Calling  routines: 
c 

c  clutr 

c 

c  Routines  called: 
c 

c  SORT 

c 

c  Glossary: 
c 

c  alpha(*)  -  ray  launch  angle  at  the  transmitter  (radians) 

c  hxmtr  -  radar  antenna  height  (m) 

c  imax  -  number  of  elements  in  alpha,  pair,  &  r  arrays 

c  psi(*)  -  grazing  angle  at  the  sea  surface  (radians) 

c  r(*)  -  range  (km) 

c  rk  effective  earth  radius  factor 

c  rmax  -  maximum  plot  range  (km) 

c 

c  Method: 
c 

c  This  subroutine  takes  radar  antenna  height  and  effective  earth 

c  radius  and  returns  launch  angle,  alpha,  and  grazing  angle,  psi, 

c  versus  range  for  values  of  psi  between  0.1  and  10  degrees 

c  using  the  solution  of  the  ray  trace  equations  in  terms  of  ae. 

c 
c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c 

SUBROUTINE  psirae  (rk,  hxmtr,  rmax,  alpha,  psi,  r,  imax) 
c 
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DIMENSION  alpha (50) ,  psi(50),  r(50) 
c 

ae  -  rk  •  8371. 

c 

c  deteraine  the  range  at  which  psi-  .1  deg 

c 

psi(l)  -  0.001745 

a  -  1. 

b  -  2.  •  ae  •  pai(l) 

c  -  -2.  •  ae  •  hxatr  /  1000. 

r(l)  -  (-b  +  SQRT(b*b  -  4.*a*c))  /  (2.*a) 

IF  <r(l)  .GT.  raax)  r(l)  -  raax 
alpha(l)  -  -(hxatr/(1000.  *r(l) )  -  r (1) / (2 .  *ae) ) 
c 

c  deteraine  range  at  which  psi  -  10  deg 

c 

b  -  2.  •  ae  •  0.1745 

rain  -  (-b  +  SQRT(b*b  -  4.*a»c))  /  (2.*a) 
c 

c  deteraine  alpha  ic  psi  for  intermediate  ranges 

c 

i  -  1 

dr  -  r(l)  -  rain 
DO  i  -  2,  45 

r(i)  -  r (1)  -  dr*i/45. 

alpha(i)  -  -(hxmtr/(1000. *r(i) )  -  r(i)/(2. *ae) ) 
pai(i)  -  hxmtr/(1000.*r(i))  -  r(i)/(2.»ae) 

END  DO 

iaax  -  45 

c 

c 

RETURN 

END 

c 

c  Subroutine  push 
c 

c  PUSH  stores  elements  in  an  array  and  returns, 
c 

c  Variable:  Description: 

c 

c  array  iq  array  to  hold  data  eleaents 

c  lq  Size  of  data  array 

c  nq  Number  of  data  elements  stored  in  data  array 

c  data  The  data  element  to  be  stored 

c 
c 

SUBROUTINE  push ( ARRAY . iq.nq, data) 
c 

real*4  data, array 

integer*2  iq.nq 
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dimension  array (iq) 

c 

c  Shift  array  elenents  down  one 

c  do  i«iq,2,-l 
DO  j-2.iq 
i-iq-( j-2) 
array ( i ) -array ( i -1 ) 

END  DO 
C 

c  Store  new  data  eleaent  in  top  of  array 

array ( 1 ) -data 
nq-nq+1 

IF(nq  .GT.  iq)  nq  -  iq 

RETURN 

END 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c  Name: 

c 

c  SUBROUTINE  upgo  (alpha,  ixatr,  idelta,  r) 

c 

c  Purpose : 
c 

c  Ray  trace  for  upgolng  raya,  alpha  >-  0 

c 

c  Inputs: 
c 

c  alpha,  idelta,  rhl(*),  rol(*),  ixatr 

c 

c  Outputs: 
c 

c  r 

c 

c  Calling  routines: 
c 

c  psir 

c 

c  Routines  called: 


c 


c 

A 

SORT 

w 

c 

A 

Glossary: 

c 

C 

alpha 

-  launch  angle  in  radians 

c 

idelta 

-  index  of  level  in  profile 

for 

evap.  duct  height 

c 

ixatr 

-  index  of  level  in  profile 

for 

transaitter 

c 

rhl(*) 

-  height  (a) 

c 

ral(*) 

-  aodified  refractivity  (M) 
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-  ranee  (ka) 


c  r 

c 

c  Method : 


c 

c  An  upgoing  ray  Is  traced  through  the  successive  layers  of 

c  the  M  profile  to  deteraine  the  range  at  which  it  reaches  the 

c  height  of  the  evaporation  duct, 

c 
c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

SUBROUTINE  upgo  (alpha,  ixatr.  ldelta,  r) 
c 

DIMENSION  rhl(10) ,  ral<10) 


C0M40N  /RAYTRACE/  rhl,  ral 
i  •  ixatr 

DO  WHILE  (i  .LT.  ldelta) 
dh  -  rhl(i  +  1)  -  rhl(i) 
da  -  raid  +  1)  -  ral(i) 
dadh-  da/dh 

ir  (dadh  .EQ.  0.)  dadh  -  le-8 

sub  -  alpha**2  +  0.000002*dadh* (rhl (i  +  l)-rhl(i)) 
IF  (sua  .LT.  0.0)  sua  -  0.0 
alphap  ■  SORT (sua) 

rp  -  r  +  (alphap-alpha)/(dadh  •  0.001) 
alpha  -  alphap 
r  -  rp 
i  -  i  +  1 
END  DO 

RETURN 

END 
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